Direct experimental test of scalar confinement 
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The concept of Lorentz scalar quark confinement has a long history and is still widely used despite 
its well-known theoretical faults. We point out here that the predictions of scalar confinement 
also conflict directly with experiment. We investigate the dependence of heavy-light meson mass 
differences on the mass of the light quark. In particular, we examine the strange and non-strange 
D mesons. We find that the predictions of scalar confinement are in considerable conflict with 
measured values. 
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I. A LITTLE HISTORY 

Within a year of the discovery of charmonium in 1974, 
a picture emerged of two massive quarks moving in the 
QCD motivated static potential pj, 



V(r) = -- 
r 



(1) 



Here the color Coulomb constant k is given by k — 4a s /3 
in lowest order perturbation theory and a is the confine- 
ment constant, or asymptotic confining force. The non- 
relativistic Schrodinger equation accounted in a natural 
way for the spin averaged cc and bb spectra and decay 
rates Q. It was soon realized that in order to under- 
stand the spin dependence of the heavy onia states one 
must further specify the Lorentz transformation proper- 
ties of the interaction. Henriques, Kellett, and Moor- 
house |3( proposed that the short range Coulombic part 
was a Lorentz vector and the long range confining part 
a Lorentz scalar. The reasoning by the above group and 
others was that in order to account for observations, 
the spin-orbit interaction must be suppressed by a par- 
tial cancelation between the short-range and long-range 
contributions. 

This picture flourished unchallenged for fifteen years 
despite the lack of success in formally relating scalar con- 
finement to QCD. Around 1990, the work of the Milan 
group p| clarified this question dramatically. Using the 
low velocity Wilson loop formalism pioneered by Eichtcn 
and Feinberg |fg, and later by Gromes [7j, they found 
both the spin-dependent and spin-independent relativis- 
tic corrections in heavy onia. Their astounding result 
was that the long-range spin-independent QCD correc- 
tions differed from those of scalar confinement, though 



the long-range spin-orbit corrections were the same pure 
Thomas ones given by scalar confinement. These results 
were subsequently verified by lattice simulations of QCD 
0. 

Because the spin-independent relativistic corrections 
to scalar confinement are incorrect, the scalar confine- 
ment scenario should logically be discarded. Indeed, the 
whole concept of potential confinement is in error. For- 
tunately, there is an alternative physical picture that can 
be employed. Back in 1982, Buchmiiller pointed out 
that a color electric flux tube should automatically yield 
the desired pure Thomas spin-orbit interaction because 
there is no color magnetic field in the quark rest frame. In 
1992, it was shown that the spin-independent relativistic 
corrections 

.10] 

to the flux tube model exactly matched 
the QCD predictions 0- Recently, we have constructed a 
consistent classical action for spinning quarks connected 
by a QCD string (flux tube) Hence there exists a 

simple physical picture that is consistent in many ways 
with QCD and is not a simple potential model. 

Although inconsistent with QCD, scalar confinement 
models remain popular, probably because of their relative 
ease of solution. In this paper we point out an instance 
where the predictions of spin-independent scalar confine- 
ment disagree directly with experiment. We consider 
here the dependence of heavy-light (HL) meson masses 
on the light quark mass. We discuss a general HL poten- 
tial model wave equation in Sec. Inland exhibit a simple 
analytic perturbative solution for scalar confinement. We 
demonstrate the high accuracy of our numerical solutions 
by comparing them to the analytic ones. In Sec. IIIII we 
collect and discuss the experimental data which we will 
compare to our predictions. In particular, we will estab- 
lish that spin splittings do not depend on the mass of the 
light quark. Comparison of three confinement scenarios 
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to data is given in Sec. II VI where we give theoretical val- 
ues obtained by exact numerical solution of the spinlcss 
Salpeter equation for both scalar and time component 
vector (TCV) confinement as well as flux tube confine- 
ment. Our conclusions in Sec. are that scalar con- 
finement predictions of the light quark mass dependence 
are not good and can be clearly seen in the difference 
of S'-wave and P-w&ve meson states. We note further 
that electric (TCV) confinement also does not account 
for this difference but that the flux tube model exhibits 
remarkable agreement with experimental data. 



II. SCALAR CONFINEMENT EQUATION: A 
PERTURB ATIVE SOLUTION AND AN EXACT 
NUMERICAL SOLUTION 

A. The heavy-light potential model 

We consider a spinless quark of mass m that moves in 
a potential field. We assume that this field consists of 
central Lorentz scalar and time component vector fields, 
S(r) and V(r) respectively. The expected Lagrangian for 
this system is 



L = -m(r)vl — v 2 — V(r) 
m(r) — m + S(r). 

The momentum and Hamiltonian arc then 

p = 771(7") 7 v, 
H = vp-L, 



where 7 = l/yl — v 2 . 

By squaring p, we see that 



m(r) 7 = \Jp 2 + m(r) 2 . 
Substitution into the Hamiltonian (J5J then yields 



H = v/p 2 + (m + %)) 2 + V(r). 



(2) 
(3) 



(4) 
(5) 



(6) 



(7) 



For a state of definite angular momentum £ the mo- 
mentum square can be written 



p 2 =Pr 2 + 2 ■ 

„2 . 



(8) 

With the usual replacement p 2 — > — \^r, we obtain 
a wave equation that can always be solved numerically 

E3. 



B. A perturbative estimate of the light quark mass 
dependence 



short range interaction, V(r) = — jjy = — From the 
general Hamiltonian 0, we find 



H = 1/p 2 + (m + ar) 2 . 



(9) 



We wish to treat the short range parameter k and the 
light quark mass m as small perturbations. In the case 
where both k and m vanish, the zeroth order Hamiltonian 
is 



Ho = VP 2 + {ar) 2 - 



(10) 



The Hamiltonian Hq has the same eigenstates as its 
square which, with the replacements of Eq. (|SJ| and p 2 — > 



1 d z 

r dr 2 



leads to the harmonic oscillator equation, 



cPuq 
dr 2 



£(£ + !) 



2 2 
— a r 



where in the limit of small r, uq 
malized to 



dr\u \ = 1. 



u = 0, (11) 
r e+1 and uq is nor- 



(12) 



The solution for the wave function and eigenvalue is stan- 
dard; 



uo(r) = 

Kt = 



Af n4 r* +1 e-i^ 



2a" 



(n - 1) 



Ei = 2a 



T{£ + n+\) ' 
1 
2 



2n- 



(13) 
(14) 

(15) 



Here L^l^(ar 2 ) is the usual Laguerre polynomial, n is a 
positive integer starting with 1, and I is a non- negative 
integer starting with zero. 

We now determine the effect of tur ning on k and m 
using the Feynman-Hellmann theorem [13j, 

dE/d\= (dH/dX), 



(16) 



where A is any parameter of the Hamiltonian. Taking 
the expectation values using the k — m — wavefunc- 
tions l)13|l will yield an expansion in these parameters. 
To leading order we have 



E = E -k{r-i) + ^(r) 



(17) 



The expectation values are worked out in general in 
the Appendix. Here we only consider the m and k de- 
pendence of the 15 and IP states (i.e., n = 1 and I = 
and 1). The results are, 



We next consider a simple model with scalar confine- 
ment, S(r) = a|r| = ar and a time component vector 



'is 



= 2«/° 



(18) 
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(19) 
(20) 
(21) 



Using Eq. 1)17(1 . we obtain the m dependence near k 
m = Q (m slope) for the energies, 



dE 



is 



dm 

d{E 1P - E 1S ) 
dm 



/3tt 

2/4 



/3?r 



0.6515, (22) 
- 1 J = 0.0214. (23) 



Again using Eq. (|17|l we obtain the k dependence near 
k — m — [k slope) for the energies, 



v/a dfc V 7r 



1.128, (24) 



1 d(.E 



ip 



si 



2 /l 



3 V 7T 



0.376. 



(25) 



C. Consistency of the exact numerical solution and 
the analytic perturbative solution 

Before we proceed to a detailed comparison of the 
predictions of scalar confinement with the experimental 
data, we pause to verify the accuracy of our numerical 
method. In particular, we compare our analytic values 
for the slopes at m = k = in Eqs. l(2"2"|) to (|23|l with 
our numerical method. This step is important because a 
general analytic solution of the spinless Salpeter equation 
is not known. Only for the remarkable case of a massless 
particle and linear scalar confinement, which is equiva- 
lent to the non-relativistic harmonic oscillator, can one 
obtain an analytic solution. In the general case, one must 
rely on exact numerical solutions. Some time ago, we in- 
troduced 0| a variational method, the Galerkin method, 
into particle physics to solve the spinless Salpeter equa- 
tion with a time component vector interaction. This very 
robust method is applicable to a wide range of differen- 
tial and integral equations. The method has been sharp- 
ened over the years by many authors |l4j . One can cope 
with eigenvalue equations for operators that are compli- 
cated functions of both momenta and coordinates, such 
as the scalar confinement Hamiltonian JJJJ, by using ba- 
sis functions that can be Fourier transformed. We have 
performed a careful numerical solution of the eigenvalue 
equation for the Hamiltonian for small k and m, and 
found the m and k slopes. The results are in excellent 
agreement with the values obtained by analytic calcula- 
tion in Eqs. through l(25)l. 



III. EXPERIMENTAL DATA 
A. Spin splitting is independent of light quark mass 

Wc first use experimental data to demonstrate, rather 
conclusively, that spin splittings within a given orbital 
angular momentum multiplet do not depend on the light 
quark mass. In particular, we consider several heavy- 
light meson spin multiplets in which both the strange 
and non-strange members have been observed |l5j . First 
we examine the D and D s type mesons with £ = (S- 
waves). The hyperfine splittings for D s and D mesons 
are 



D* s 
D* ± - 
n* 



D± 
-D 



143.8 ± 0.4 MeV, 
140.64 ± 0.10 MeV, 
142.12 ±0.07 MeV. 



(26) 



The corresponding S- wave hyperfine splittings for B type 
states are 

B* - B s = 47.0 ± 2.6 MeV, 
B* - B = 45.78 ± 0.35 MeV. (27) 

It is clear that the substitution of a strange quark for 
a non-strange light quark changes the 5-wave hyperfine 
differences by at most a few MeV. 

We next consider some measured P-wave heavy-light 
spin splittings. From |l5j we find, 



D s2 - D sl = 37.0 ± 1.6 MeV, 
Dl-Dl = 36.7 ± 2.7 MeV, 
D 2 - D i = 32 ± 6 MeV. 



(28) 



Again we note the apparent vanishing of light quark 
mass dependence this time in a P-wave spin splitting. We 
conclude that spin splittings are only weakly dependent 
on the light quark mass. We will exploit this fact in 
section Till CI 



B. The light quark mass dependence of a IP- 
difference 



IS 



In the preceding subsection we observed from exper- 
iment that both iS-wave and P-wave heavy-light spin 
splittings (within a spin multiplet) were independent of 
the light quark mass. We next consider the mass split- 
tings between pairs of states corresponding to different 
orbital angular momenta and examine the light quark 
mass dependence of this difference. We choose the D\ 
P-wave state and the pseudoscalar D meson. The best 
measurements are 



D" 



557.5 ± 2 MeV. 



(29) 



When the u light quark is replaced by a strange quark, 
the corresponding difference becomes 



D sl - D s 



567.3 ± 0.4 MeV. 



(30) 
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The differences A u and A s are amazingly similar. We 
see that they differ by 

A = A s - A u = 9.8 ± 2 MeV. (31) 

C. Conclusion for the spin-averaged IP — IS 
difference 

We demonstrated in Sec. If If Al that both 5- and P- 
wave spin splittings are, within error and isospin uncer- 
tainty, independent of light quark mass. We may there- 
fore conclude that the m-dependence of the difference 
A s — A u also represents the m-dependence of the spin- 
averaged IP— 15 excitation energies. To show this ex- 
plicitly, we write the HL meson mass as the sum of the 
heavy quark mass and the excitation energy, 

M = m Q + E. (32) 

We then separate the excitation energy into spin- 
averaged and spin-dependent parts, 

E = E SA + E SD . (33) 

As in Eqs. (|2*9"j> and (3D) , we define A s and A„ to 
be the differences between P-wave (J p = 1 + ) and 5- 
wave (J p = CP) states, for strange and non-strange 
light quarks respectively. Using our observation that the 
spin-dependent parts are essentially independent of light 
quark mass, we see that the spin dependent parts will 
cancel in the difference A = A s — A u . Thus 

A = (Eip - E 1S ) S - {E XP - E 1S ) U (34) 

measures the light quark mass dependence of the spin- 
independent IP— 15 level difference. 

Two main results emerge from this subsection. First, 
the light quark mass dependence is only about 2% of 
the measured level difference. This is a useful tool in 
the analysis of heavy-light spectroscopy 0. Second, 
although it is small, A is definitely not zero. We will use 
the actual value l|3l"|) to test the predictions of various 
assumptions about the nature of quark confinement. 

IV. COMPARISON WITH EXPERIMENT 

We are now prepared to compare carefully the different 
confinement pre dictions with experiment. As we orig- 
inally noted [0| , the m dependences of all the heavy- 
light states are amazingly similar. We have also noted in 
Sec. Ill Cl that this universal m dependence is nearly sat- 
isfied in an analytical calculation. In the example with 
scalar confinement, the m slope for the difference IP— 15 
is about 30 times smaller than each separate slope. Fur- 
thermore, we noted in Sec . 1111 Bl that when one compares 
two heavy-light states, one P-wave and one 5-wave, the 
difference changes by less than 2% when a non-strange 



light quark is replaced by a strange one. This change is 
not zero however, but for the D\ and D states has the 
experimental value (|31|) 

A = A s — A M = 9.8 ± 2 MeV. (35) 

In Fig. ^ we show the numerical mass splittings of heavy- 
light mesons for three different confinement scenarios, all 
with the same short range energy — fc/r, and the same 
k = 0.5, as a function of light quark mass. Each of 
the three confinement scenarios has the same asymptotic 
confinement force, a = 0.18 GeV 2 . The upper curve 
assumes linear scalar confinement, the middle curve is 
the prediction of the relativistic flux tube and the lower 
curve linear time component vector confinement. In the 
scalar and time component vector potentials, the poten- 
tials S(r) and the long-range part of V(r) respectively 
are ar. In the case relativistic flux tube model the string 
tension is a. 




_ 10 i , , , , , 1 

100 200 300 400 500 600 



m [MeV] 

FIG. 1: Dependence of the IP — IS energy level difference on 
the light quark mass m. The numerical results compare dif- 
ferent confinement mechanisms (all with the same long range 
confinement force a = 0.18 GeV 2 and all three calculations 
assume the same short range Coulombic constant k = 0.5.) 

We note that the scalar confinement has the most rapid 
increase of the IP— 15 difference as the light quark mass 
increases. The flux tube confinement is intermediate 
and time component vector linear confinement actually 
decreases slightly. Using the reasonable values for the 
strange quark mass (500 MeV) and the non-strange mass 
(300 MeV) [H, we find the following values for A, 

^scalar 

19 MeV, (36) 

^flux tube 

= 10 MeV, (37) 
Atcv = -1 MeV. (38) 

In comparing these values to the experimental value given 
in Eq. (ED, 

A oxp = 9.8 ± 2 MeV, (39) 
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we observe that both scalar confinement and time com- 
ponent vector confinement are quite inconsistent with the 
experimental result. However, the flux tube confinement 
prediction is in agreement with the experimental result. 
These results are depicted in Fig. |2] 



20+- 



FIG. 2: The difference in the IP— IS energy gap between 
a light quark mass value of 500 MeV and 300 MeV. The 
experimental value is the dot below the 10 MeV level. All 
theoretical predictions have the same long range confinement 
force, a — 0.18 GeV 2 , and the same short range Coulombic 
constant, k = 0.5. 




FIG. 3: The IP — 15 level difference as a function of the 
Coulomb constant k for scalar confinement. In this case there 
is a noticeable dependence on the light quark mass m and we 
show the calculation for m — 0, 300, and 600 MeV. Since the 
experimental value was found for D s mesons we see a value 
of k — 0.53 is indicated with m = 500 MeV. 



strange light quark masses. This relation is the result of 
the measured value for the B s and B mass difference, 



V. DISCUSSION 

A. m slopes and the k value 

The reader may have noticed that the m slope for the 
IP— 15 difference calculated analytically (|23[) is about 
0.02 whereas the m slope in Fig. Q] is about five times 
larger. This is because the Coulomb constant A; is zero in 
the analytical result and has the more realistic value 0.5 
in the numerical calculations shown in Fig. 2] The effect 
is magnified due to the strong cancelation in the IP— 15 
difference. 

The choice k — 0.5 is quite reasonable. As we noted in 
[l6|. the spin-averaged experimental value for E\p — E\s 
can be extracted from the D s mesons; 

E XP -E 1S = 439 MeV. (40) 

The appropriate choice of k is determined by computing 
the above quantity as a function of k as shown in Fig. 
We conclude that the common value of k = 0.5 assumed 
in the considerations of the last section was appropriate 
and that the effect of a slightly larger k would only make 
our conclusions stronger. 



B. Constituent quark masses 

Another interesting result of a choice of confinement 
scenario is a constraint that relates the non-strange and 



B 8 — B ~ 91 MeV. (41) 



600.0 




100.0 

o.o I 1 1 1 1 1 1 1 

50 100 150 200 250 300 350 400 

m u [MeV] 

FIG. 4: The constraint between m s and m u light quark 
masses due to the measured B 3 — B mass difference 1411 . 
The solid curve assumes scalar confinement while the dashed 
curve is due to flux-tube confinement. The "data" points are 
quark masses analyses of hypero n magnetic moments in the 
constituent quark model pH Il7|l . 

In Fig. 01 we show this constraint scalar confinement 
and flux tube confinement scenarios. Although the curves 
differ considerably for small light quark mass, they are 
quite similar in the larger range. On the figure we show 
constituent quark masses obtained by analyses of hy- 
peron magnetic moments |l5l llT^ . The result justifies 
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our choice of 500 MeV for the strange quark and 300 
MeV for the non-strange quark. 



VI. CONCLUSIONS 

Our points in this paper may be quickly summarized. 

• From the data, we conclude that spin splittings do 
not vary with light quark mass value. We then can 
extract the change (|3 1|) in the spin-averaged IP— IS 
level difference as one replaces the non-strange by 
a strange light quark. 

• We compare the predictions of the scalar poten- 
tial, time component vector potential and flux-tube 
quark confinement scenarios with experimental re- 
sults, as shown in Fig. [21 The conclusion is that 
flux-tube confinement works well while both scalar 
and time component vector confinement fail badly. 

We observe therefore that scalar confinement has at 
least one specific point of disagreement with experiment. 
This complements the theoretical disagreements with 
QCD mentioned in the introduction. 
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APPENDIX A: EXPECTATION VALUES 

The principal aim here is to compute the expectation 
value of r p with the harmonic oscillator wavefunctions 
(E3), 



With the choices 



poo 

r p )=Nl e / drr 2e+2 +Pe- ar2 
Jo 



ifcOxr*) . (Al) 



A change of integration variable to the dimensionless 
combination z — ar 2 yields 



a = £+-, 



(A6) 
(A7) 
(A8) 



we substitute Eq. (|A3jl into \K2\ and use the orthonor- 
mality relation for Laguerre polynomials |Tsj 



dz z? e~ z Lj (z) 4 (z) = F(j+ ,f + 1) 5„, t A!) ) 



3- 



to obtain our general result, 

-, 2 



a 



E 



P 

2/ n—3 



n-i\ nj + e + k + i) 



(A10) 

In Eq. ljAT0)l we use the notation for the binomial coeffi- 
cients, 



(AH) 



m J m\(n — m)\ 



The specific result that is required in Sec. IIIBl is for the 
ground state (n = 1) is 



\ / n—1 



Some cases of direct interest arc 



p=-l 
p = 0: 
p=l: 

p = 2: 



(r- 1 ) = 
(1)=1, 



. r(i + i) 
r(^ + |)' 



i r(£ + 2) 



(A12) 

(A13) 
(A14) 
(A15) 

(A16) 



( r P)=N 2 



L n-\( Z ) 



(A2) 

It is helpful to use the Chu-Vandcrmonde sum formula 

M, 



^% Lf_xW, (A3) 



where the Pochhammer symbol (z)n is defined as, 

(z) N = z(z + l)---(z + N-l) = E(f+p,(A4) 

(z) = 1. (A5) 



[1] E. Eichten, K. Gottfried, T. Kinoshita, J. Kogut, K.D. 
Lane, and T.-M. Yan, Phys. Rev. Lett. 34, 369 (1975). 

[2] Quarkonia (W. Buchmiiller, ed.) Current Physics- 
Sources and Comments Vol. 9, (North-Holland, 1992). 

[3] A.B. Henriques, B.H. Kellet, and R.G. Moorhouse, Phys. 
Lett. 64B, 85 (1976). 

[4] N. Isgur, Phys. Rev. D 57, 4041 (1998). 

[5] N. Brambilla and G. M. Prosperi, Phys. Lett. 236B, 69 
(1990); Phys. Rev. D 46, 1096 (1992); 

[6] E. Eichten and F. Feinberg, Phys. Rev. D 23, 2724 
(1981). 

[7] D. Gromes, Z. Phys. C 26, 401 (1984). 
[8] G. Bali, Phys. Rept. 343, 1 (2001). 
[9] W. Buchmiiller, Phys. Lett. 112B, 479 (1982). 
[10] M.G. Olsson, C. Olson, and Ken Williams, Phys. Rev. D 

45, 4307 (1992). 
[11] T.J. Allen, M.G. Olsson, and J.R. Schmidt, Phys. Rev. 

D 69, 054013 (2004). 
[12] Steve Jacobs, M.G. Olsson, and Casimir Suchyta III, 



Phys. Rev. D 33, 3338 (1986); M.G. Olsson, S. Veseli 
and K. Williams, Phys. Rev. D 51, 5079 (1995). 

[13] H. Hellmann, Einfuhrung in die Quantenchemie 
(Deuticke Verlag, Leipzig, 1937); R.P. Feynman, Phys. 
Rev. 56, 340 (1939). 

[14] L.P. Fulcher, Phys. Rev. D 50, 447 (1994); C. Semay 
and B. Silvestre-Brac, Phys. Rev. D 52, 6553 (1995); 
W. Lucha, F.F. Schoberl, Int. J. Mod. Phys. A 15 3221 
(2000); Phys. Rev. A 60, 5091 (1999). 

[15] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002) and 
2003 off-year partial update for the 2004 edition available 
on the PDG WWW pages (URL: lhttp://pdg.lbl.gov/t 

[16] T.J. Allen, Todd Coleman, M.G. Olsson, and Sinisa 
Veseli, Phys. Rev. D 69, 074010 (2004). 

[17] Jerrold Franklin, Phys. Rev. D 66 033010 (2002). 

[18] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Se- 
ries, and Products (Academic Press, New York, 1980). 



